Quantum Zeno stabilization in weak continuous measurement of two qubits 
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We have studied quantum coherent oscillations of two qubits under continuous measurement 
by a symmetrically coupled mesoscopic detector. The analysis is based on a Bayesian formalism 
that is applicable to individual quantum systems. Measurement continuously collapses the two-qubit 
system to one of the sub-spaces of the Bell basis. For a detector with linear response this corresponds 
to measurement of the total spin of the qubits. In the other extreme of purely quadratic response 
the operator a^cTy + o~\o% is measured. In both cases, collapse naturally leads to spontaneous 
entanglement which can be identified by measurement of the power spectrum and/or the average 
current of the detector. Asymmetry between the two qubits results in evolution between the different 
measurement subspaces. However, when the qubits are even weakly coupled to the detector, a kind 
of quantum Zeno effect cancels the gradual evolution and replaces it with rare, abrupt switching 
events. We obtain the asymptotic switching rates for these events and confirm them with numerical 
simulations. We show how such switching affects the observable power spectrum on different time 
scales. 

PACS numbers: 73.23.-b; 03.65.Ta; 03.65.Xp; 03.67.Lx 
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I. INTRODUCTION 

The quantum Zeno effect is an intriguing prediction 
that arises in the context of quantum measurement 
theoryji* 2 ^ It states that a sequence of strong orthodox 
measurements can "freeze" a system in its quantum state, 
so that in the Zeno limit of very frequent measurements, 
the system is prevented from decay and/or subsequent 
evolution. Since its original formulation^^ there has 
been continued theoretical developmentSi2iS as well as 
experimental demonstrationist in an ensemble of quan- 
tum micro systems - trapped ions (or atoms) probed via 
fast (strong) quantum-optical measurements. 

With the recent experimental advances in the delib- 
erate fabrication of two- level quantum systems, such 
as superconducting qubitsi^ or double quantum dotsi^ 
(DQD), it has become possible to perform experiments on 
individual quantum systems (see, e.g., Refs. 26«lfi). How- 
ever, it can be very difficult to make repeated strong or- 
thodox measurements on such systems. Instead, a more 
practical measurement scheme utilizes a detector such 
as a quantum point contact (QPC) or a single-electron 
transistor 17 (SET) that is weakly coupled to the qubit(s). 
The measurement record in such a situation is a fluc- 
tuating current I(t) that accumulates a distinguishable 
signal-to-noise ratio after some timeiZ^ 

On one hand, the weak coupling between detector and 
qubit(s) permits the quantum system to remain relatively 
well isolated from "outside" classical noise. On the other 
hand, it means that instead of simple abrupt collapse, 19 
we have to deal with a theory of continuous (weak) mea- 
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surements of a single quantum system. The development 
of such theories started long ago~^i^ and has attracted 
most attention in quantum opticsi^Si 2 ^ Despite the sim- 
ilarity of their underlying principles, the theories may dif- 
fer significantly in formalism and area of application; for 
solid-state qubits a "Bayesian" theory was developed rel- 
atively recently^ This formalism is reviewed in Ref. EtI 
and its equivalence to the quantum trajectory approach 
translated^SiS 9 . from quantum optics is shown in Ref. l28l 

For single-qubit continuous measurement ^22*^1 it was 
shown that in a hypothetical situation of very strong con- 
tinuous qubit-detector coupling, it would be possible to 
approach the quantum Zeno limit. For large (but finite) 
qubit-detector coupling long periods of "freezing" of the 
single-qubit state are interrupted by rare, abrupt jumps 
between qubit states. Such jumps produce corresponding 
jumps in the average detector current signal that closely 
follows the qubit evolution. This is in agreement with 
the classical character of the strong coupling measure- 
ment regime jSLiSS 

In this paper we show that the Zeno-like regime can 
arise even in a relatively practical weak continuous mea- 
surement context. We initially consider two identical 
qubits coupled symmetrically to a detector. Measure- 
ment continuously collapses the two-qubit state into mea- 
surement subspaces, i.e., stable "points" or stable multi- 
dimensional subspaces of the two-qubit Hilbcrt space. 
On a time scale greater than the characteristic measure- 
ment time the collapse is like a projective measurement 
in the Bell basis. For a detector with linear or non-linear 
response^ this corresponds to measurement of the total 
"spin" of the qubits (if^ + l?^) 2 . In the degenerate 
case of purely quadratic detection^ the spin-1 subspace 
itself splits and effectively the operator <TyO-y' +cr^a^ 
is measured. The Zeno-like regime manifests itself in 
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a situation of slightly asymmetric qubits, wherein the 
Hamiltonian attempts to evolve a collapsed state out 
of its measurement subspace. Our results show that 
measurement stabilizes the state, stifling this evolution. 
Instead of gradual evolution, one finds rare quantum 
switching events between different subspaces. 

Given the noisy experimental signal I(t), such switch- 
ing events can be identified by consideration of the power 
spectrum and/or average current. Using the stochastic 
Bayesian approach we evaluate the power spectrum in 
several important cases including the linear and purely 
quadratic detection regimes. We estimate the measure- 
ment time needed to accumulate the spectrum and show 
that even for a non-ideal detector it is possible to observe 
switching events if the lifetime of the subspaces is large 
enough. On a long time scale compared to the subspace's 
lifetime, the switching events lead to an averaging of 
the power spectrum (effective decreasing of the spectral 
peaks) and also may produce a telegraph noise peak near 
zero frequency. The resulting long-time-average power 
spectrum is in accordance with both numerical simula- 
tion of the spectrum through Bayesian equations^ and 
with the ensemble averaged approachi25i2£i 

As expected in the quantum Zeno effect, the state "sur- 
vival probability" evaluated from quantum master equa- 
tions decays quadratically in time 4,5 ' 6 for small times and 
exponentially for timesi that are sufficiently large com- 
pared to the measurement time. We have obtained an- 
alytically the corresponding "lifetimes" of the subspaces 
for small deviations from the symmetric situation. We 
have verified our results by Monte Carlo simulations of 
the measurement process. 

The presentation of our results is organized as follows. 
In Sec. II, we review the Hamiltonian and the Bayesian 
equations of motion for weak continuous measurements of 
a two-qubit system. Sec. Ill describes in detail the mea- 
surement collapse scenarios that result from the Bayesian 
equations. In Sec. IV, the quantum Zeno effect is demon- 
strated when the two qubits are asymmetrical. Sec. V 
presents conclusions. 



II. MODEL HAMILTONIAN AND BAYESIAN 
EQUATIONS OF MOTION 

A. Model Hamiltonian 

Fig. nja) shows possible realization of the setup that 
we analyze. Each qubit is made of a double quantum 
doli 1 ^ (DQD), occupied by a single electron, while the 
detector is a quantum point contact (QPC) located in be- 
tween DQDs. Another possible realization (not shown) is 
based on single-Cooper-pair-box (SCPB) qubits— mea- 
sured by a single-electron transistor (SET). 

In the Hamiltonian of the system, 

Ti- = Hqb +'Hdet + ~Hint, (f) 




FIG. 1: Schematic of two qubits measured by a linear and 
non-linear detector, (a): Realization based on double quan- 
tum dots measured by a QPC. (b): Possible realization of 
a purely quadratic detector based on double quantum dots 
measured by two QPCs. 

the first term describes two qubits alone with direct in- 
teraction between them, 

Kqb = (e Q /2)(a|a i - aja T ) + H a (a\ai + a{a T ) 
+ (e b /2)(b\b l - 6|6 T ) + H b (b\b L + 6{6 T ) 

+U(a\ai-a\a l )(b\bi-b\b l ) (2) 

The amplitudes H a and describe tunneling within 
each qubit while e a and Eb are the energy biases of each 
qubit. In the realization shown in Fig. ^a), U gives the 
strength of the Coulomb interaction. For other realiza- 
tions, the interaction may arise from a different physical 
source. In the absence of interaction (U — 0), each qubit 
becomes an independent two level system. The free Rabi 
oscillation frequencies of qubits "a" and "b" are then 
given by h a = (4H 2 + e 2 a ) 1/2 /h and h b = (AH 2 + e 2 )^ 2 /ft 
respectively. In the following it will be convenient to de- 
fine the Rabi frequencies at zero energy bias £l a = 2H a /h 
and SI;, = 2Hb/h. 

As an example of a detector we consider the case^I of 
a low transparency QPC, with Hamiltonian 

HdeT = E l°h +J2 E r4cr + J2( Tc l Cl +Hx -) ( 3 ) 

I r (,r 

An operator of the form c\ (c£) creates an electron in the 
left (right) lead of the detector. The final term causes 
electrons to make transitions between the leads. We sum 
over all states I in the left lead and all states r in the 
right lead. The qubit-detector interaction term can be 
written as 

Hint = ^(aja T - aja|)(6j6 T - &}&;)x4c; + He. 

l,r 

+ y](o}of - ajaj.) AT a c£q + H.c. 

l,r 

+ 5>j&T - AT b c tci + H.c. (4) 

l,r 

The amplitudes T and AT ai ATf,, x are assumed inde- 
pendent of the energy of the tunneling electrons. For 
simplicity, we will consider a situation when there is no 
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relative phase between these amplitudes. Such a relative 
phase would lead to additional decoherence i—————^ From 
the form of (@J, one sees that the detector's measurement 
basis is jf), \\) for each qubit. Eq. Q is the most gen- 
eral detector-system interaction in this basis in the case 
of two qubits. 35 Higher terms in Eq. (@J will appear for 
larger number of qubits. 

Purely quadratic detectors (x / and AT a ^ = 0) 
may be of use to perform simple quantum error correc- 
tion protocols— To realize such a detector physically, 
one possibility is to form two identical QPCs into two 
arms of an Aharonov-Bohm type loop (no magnetic field 
is applied). If two qubits are placed exactly at the geo- 
metrically symmetric point of the loop as in Fig. IHb), we 
have a purely quadratic detector which responds only on 
the relative state of the two qubits, i.e. it does not distin- 
guish between states |totfe) an d ||alf>) and also does not 
distinguish between |fa4&) and |j. T&). Ref . l36l presented 
an explicit example of a non-linear detector constructed 
to measure flux qubits in the form of a superconducting 
loop interrupted by two Josephson junctions. 

The four basis states of two qubits, |1) = |TaT&), |2) = 
ItaJ-b), 1 3) = |4aTb>, |4) = \iaib), correspond to four val- 
ues of the average current through the detector: 

ii, 2 = 2vr(T ± X + AT a ± AT b ) 2 P i Pr e 2 V/h 

I 3A = 2vr(T T X - AT Q ± AT b ) 2 PlPr e 2 V/h, (5) 

where V is the QPC voltage and pi and p r are the den- 
sities of states in the left and right leads. 

Note that if a very large voltage were applied to the 
QPC for a short time, the current variances would be- 
come much smaller than the corresponding current dif- 
ferences. This is probably experimentally impractical, 
but in principle it would constitute a strong orthodox 
measurement. 



B. Weak continuous measurements of qubits 

In the course of a measurement, our knowledge about 
a quantum system is updated according to the result 
of the measurement. In the case of strong instan- 
taneous measurements the update ("collapse") is also 
instantaneous^— In a weak continuous measurement, the 
change of the system's quantum state occurs gradually. 26 
The measured result is intrinsically noisy (it is a noisy 
quasi-continuous current I(t) or voltage V(t) in the sim- 
plest solid-state realizations), so that the information 
from the detector is also acquired only gradually in time. 

Quantitatively, such a random measurement process 
can be described by Bayesian equations—?- for the system 
density matrix pij (t) that take into account the actual 
measurement result I(t). For a system of qubits mea- 
sured by a common detector they can be written!— 1—— as 
stochastic equations (in Ito form— £): 



(W,) 2 
45 



Hi 



Pij - T[H-QB,p]ij, 



(6) 



where the summation is performed over all of the states 
of the system, the Ik are the corresponding current levels, 
So gives the background low- frequency noise in the detec- 
tor, and 7y is the decoherence rate. For a 2-qubit system, 
there are 2 2 = 4 states in the sum over k and the Ik are 
given by JSJ). Eq. ijfjj assumes that the detector is func- 
tioning in a weakly responding regime: |Zj — I A -C Ii, Ij', 
this justifies the use of a single low frequency noise Sq for 
all states |A;) of the system. 

Eq. © can be obtained by a general "informational" 
Bayesian approach— ————i (see also an analogous deriva- 
tion of the Bayesian equations for the measurement of a 
nano-resonatofl—L). Since QPC is an ideal detector with 
100% quantum efficiency and no "phase asymmetry" , the 
evolution due to measurement basically reduces to the 
"quantum Bayes theorem."— It says that at a sufficiently 
short time scale in which the evolution by TLqb can be 
neglected: 1) evolution of the diagonal density matrix el- 
ements in the measurement basis coincides with evolution 
of classical probabilities described by usual Bayes formula 
and 2) the purity measure pij / y/pupjj stays constant. 
Such a derivation does not assume linearity of the mea- 
surement device, so it is also relevant for the general non- 
linear coupling of Q. Essentially, Eq. JBJ describes the 
interplay between the "intrinsic" qubit evolution (Rabi 
oscillations) and the evolution due to measurement. It 
implicitly assumes that the internal detector dynamics 
of |JSJ and Q is much faster than the qubit dynamics 
due to Hqb and measurement. This is seen explicitly 
in another approac h 26 ''?™ that leads to the same Eq. . 
It uses generalized Bloch equations that couple the qubit 
dynamics and the number n of electrons that have passed 
through the detector (for the one-qubit case see Ref.©. 
One assumes that the typical detector decoherence time 
is much smaller than the characteristic electron tunneling 
time, U/eV <C e/Jo so that coherences between different 
passages of electrons can be neglectedi— ^ In order for the 
current to be considered continuous, the characteristic 
electron tunneling time must be much less than the time 
scale of the qubit dynamics, e/io <C 2tt/H. Then, in 
the weakly responding regime, when individual tunnel- 
ing events produce only a small change of the qubit state 
it is convenient to condition the evolution on the current 
I(t). These conditions are sufficient for the stochastic 
process Pijit) to be Markovian and make it possible to 
obtain equation 

To simulate individual realizations of the random mea- 
surement process, the noisy detector current I(t) can be 
expressed as 



/(t) = £(i) + 5>fc fc 4, 



(7) 



Pij 



[i{t) - pM h + ij 



k J 



Pil 

So 



where £(t) represents white noise with spectral density 
So- For a measurement with an ideal detector there are 
no individual qubit decoherences (7^ =0). As a re- 
sult, Eq. I© leads to purification of the qubit state in the 
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course of continuous measurement even though the qubit 
undergoes random time evolution^ 

Averaging over different realizations of the measure- 
ment process I(t) in Eq. (jSJ means total ignorance of the 
particular measurement result and is equivalent to en- 
semble averaging. The ensemble average of the qubit 
state does exhibit decoherence even for an ideal de- 
tector. Averaging in Eq. 10 is simply achieved by re- 
moving the noise termjH The result is the standard 
ensemble-averaged equations 4 ^ with ensemble dephasing 
rates T lj = {U - Ij) 2 /4So + 7y 



Linear response 

— /(Tt) 



M 



Nonlinear response 
/(tt) 



Quadratic response 



AI+SI 



■ /(ti)=/(4t) 



/(U) 



/(ti)=/(4t) 
/(U) 



/(U)=/(lt) 
/(U)=/(tt) 



FIG. 2: Current levels for symmetrically coupled qubits, 
AI a = A/;, = AI. Linear, non-linear, and purely quadratic 
response. 



Pij — Fy pij 



-r[HQB,p]ij- 



(8) 



We see that even in the case 7y = the averaging has 
produced a decoherence rate (7, — Ij) 2 /4Sq which is the 
quantum limited decoherence from an "ensemble aver- 
aged" point of view^i 16 ! 46 ! 47 

Individual qubit decoherence 7y 7^ will arise if the 
detector is not ideal and may be understood as partial ig- 
norance of the measurement result. 26 The contribution to 
"fij is parameterized by the detector ideality or efficiency 
77 (0 < 7] < 1) 



lij = 



(9) 



The partial ignorance may arise due to extra output 
noise 4 ^ or when the QPC detector has a finite tempera- 
ture TiS& For finite temperature, one takes into account 
the total current through the QPC, I(t) = J— («)-!,_(*), 
but ignores information contained in the partial currents 
"to the right" and "to the left" in the detector. Corre- 
spondingly the quantum ideality (efficiency) is reduced 
down to^a t) = tanh 2 (eV/2T)M Notice that no mu- 
tual dephasing need be present between two states i,j if 
/, = Ij. Indeed, in this case the detector cannot distin- 
guish between the two states and therefore there is no 
(relative) information to be received and subsequently 
ignored. Then both 7^ and Tij vanish. Of course, ad- 
ditional qubit decoherence can arise due to individual 
coupling of qubits to the environment. The equation for 
7ij will then change as we will see in Sec. IV. 



III. SYMMETRIC WEAK CONTINUOUS 
MEASUREMENTS OF TWO IDENTICAL 
QUBITS 

In the quantum Zeno effect, measurement suppresses 
Hamiltonian evolution, repeatedly collapsing the system 
back to a given state. To understand how this suppres- 
sion happens, it is essential to develop an understanding 
of the process of collapse. In this section, we consider 
symmetric weak continuous measurement of two identical 
qubits. In this case the system collapses into measure- 
ment subspaces that can be identified and characterized. 
We present a detailed analysis of the possible collapse 



scenarios for coupling to a linear, non-linear, or purely 
quadratic detection device. This sets the stage for the 
appearance of quantum Zeno physics in Sec. IIVI 



A. Preliminary view of dynamics 

To consider solutions of it is illuminating to an- 
alyze first the dynamics of the Hamiltonian J5J alone. 
Introducing the Bell basis proves convenient: 



\i) B = (|u 6 >- U„T6»/V2, 
|2) s = (|T a T fc )-|U b »A/2, 

|3) B = (\UU)+ \UW)/V2, 
\4) B ^(\Uh)+\U b ))/V2. 



(10) 

(11) 

(12) 
(13) 



Using Eq. J2J it is straightforward to obtain the free evo- 
lution of the two qubits when the qubits are unbiased, 
£ a = £6 = 0, and do not interact, U = 0. In this case, 
the "+" and "— " subspaces of the Bell basis do not mix: 

|1) B -» coe(Aftt/2) |l) s -iwx(Mlt/2) \2) B (14) 
|3) B -> cos(fii) |3) B - i sin(fif) |4) B , (15) 

where AJ7 = f2 a — fib and O = (f2 a +f2j,)/2. This property 
persists when the direct interaction is switched on 



|l) s [cos(n-t) + t-^p sin(Q-i)] \1) B (16) 



u 



(17) 



|3) 



B 



u 



[cos(fl+t) + i-— sin(fi + i)] |3) 



B 



n 



-1— sin(0+i)|4) J 

n~ = V(A^/2) 2 + (u/h) 2 
n+ = V^ 2 + (u/ny. 



(18) 

(19) 
(20) 



Notice that |l) s and \2} B are eigenstates oITLqe if AO = 
il a — flf, = 0. When the energy biases e a , Eb are non-zero, 
there is no simple solution to the free equations of motion. 

Turning to the coupling Hamiltonian , we note that 
our assumption of symmetric qubit couplings implies 
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AT a — AT b . In this case, two currents in (J5J coin- 
cide, I 2 — I3 = I23, as seen in Fig.|3 so the measure- 
ment cannot distinguish between the states |2)=|t a j;,) 
and |3) = |! a t&)- Qualitatively, the measurement alone 
tends to continuously collapse the system towards one 
of the subspaces with definite current. For a linear de- 
tector (x = in Eq. |0J), there are three possible values 
of current I±, / 4 , 123. Therefore, measurement collapses 
the system toward |1) or |4) or the subspace {|2), |3)}. 
The current values are evenly spaced I\ — I23 = I23 — I 4- 
For a non-linear detector (\ 7^ in 10}), there are again 
three possible currents, so that measurement again col- 
lapses the system toward the same three subspaces. How- 
ever, the currents are not evenly spaced. For a purely 
quadratic detector, (AT atb = in |@J while \ 0), 
Ii = I4 and I2 =13, so the measurement itself collapses 
the system toward the subspace {|1), |4)} or toward the 
subspace {|2), |3)}. 

To describe the measurement process in more detail, it 
is convenient to introduce dimensionless qubit-detector 
couplings C a = h(AI a ) 2 /S H a and C b = k{AI b ) 2 /S H b 
where the current differences are defined as AI a = I2—I4, 
AI b = I3 - I4, and AI = (AI a + AI b )/2; weak cou- 
pling corresponds to C a ,C b < 1. To characterize the non- 
linearity of the detector, it is useful to define 



5I_ = [(h - h) - (I 2 - h)} 
AI [(J 2 + 1 3 )/2 - h] ' 



(21) 



/0 



V 



P12 



P23 



Pl2 + P?i \ 

P43 

P32 + Pu 







The entries indicated with a dot are determined by the 
hermiticity of the matrices. 
The noise term itself reads 



^j-term = — £(t) x 



S 

ZPll z Pl2 + 2 
_L 
2 



(24) 



Z(P22 ~ \) Z P23 



Pi?, 
2 



ZP24 



/o 



9X2 



P?4 



P22 ~ P44 P23 + P43 2 P: 



7n B B _ P12 

z Pl3 z Pl4 2 

B , B 
P22 ' P44 

2„ 

B B P23 

Z P33 Z PZ4 ~ 2 

Z (P44 -3) ) 

PX2 + Pu \ 
24 + (P22 + P44) 
Pi 2 + P34 







V • - (P22 - P44) J 

where it proves convenient to introduce the quantity 

A f ,„ 
2/ 2 



z = pu (1 + S) - p 4 4 = 2Repf 4 1 



{P22 



Pi\) 



(25) 

given in the measurement basis and in the Bell basis. 
Since the detector signal Q can be expressed, using l(2*T|) 
and <(23|l as, 



where SI = — (I2 — Ii)- These definitions apply in 

general to the non-symmetric coupling case; for identical 
symmetrically coupled qubits, AI a = AI b = AI and 
C a = C b = C, and SI is indicated on Figd For linear 
coupling (x = 0) one has 8 — while for the purely 
quadratic case 5 = — 2. 



B. Bayesian equations in the Bell basis 

It is instructive to rewrite the Bayesian equations © 
in the Bell basis flU} - (El 



I(t) 



'23 



■ AI z(t) + 



(26) 



we see that z(t) gives the two-qubit signal. In addi- 
tion, for linear coupling z(t) is equal to the average z- 
component of the total spin. 

The unitary part has one contribution proportional to 
the Rabi frequency f2, one proportional to the energy bias 
e, and one proportional to the interqubit coupling U /h: 



dpf 3 
dt 



H- 



dpf 3 
dt 



, d 4 


H 




dt 

a 


£ 


dt 



(27) 



where 



±1 
dt 



dpf 3 
dt 



dpi 
dt 



(22) 



H- 



where we distinguish a non-unitary part due to measure- 
ment from a unitary part due to free Hamiltonian evolu- 
tion. For qubits with identical couplings, these contribu- 
tions (in Ito form) are: 



dp% 
dt 



/0 



£ij-term 



Pl2 



(A/) 2 




(23) 



Pl4 



^m(2 + «5) 2 pg iImpf 4 (2 + 5) 2 







P34 



P44-P22 ( 2 + S) 2 J 



dpi 
dt 



dfl 
dt 



and 



= O 



/0 
• 

v- • 



/o 



dpi 
dt 



ipf 4 

iP24 



iPi3 

¥23 



-2Impf 4 i(pg- P g 
2Impf 4 



/ 



ip?4 
-2Impf 4 



¥12 

°22 - P44) 
I) ipg 

'24 





-iP43 AP2 
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¥12 



zpf 4 \ 

-ipg 

ipf 4 
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(28) 



2Impf 4 J 



(29) 



(30) 
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The state (f 101) which we have denoted has spin-0, 
while the other three Bell states (|llf) - <|13ll comprise the 
spin-1 subspace. By inspecting Eqs. and J2U, one 
sees that measurement does not mix the spin-0 and spin- 
1 subspaces. Indeed, when the system is in the spin-0 
subspace, pf x = 1, then dpfi/dt vanishes (since z = 0). 
It also vanishes when the system is in the spin-1 subspace 
(pfi = 0). Since Trp B = 1, it follows that if the system 
is initially in the spin-1 subspace it will remain there. 
In the purely quadratic measurement case (5 = —2) we 
have from 1(2*51) z = — (p B 2 + p B 4 ), and calculation shows 
that the spin-0 state |l) s , the state |2) B , and the "+" 
subspace, {|3) B , |4) B }, similarly do not mix. For unbi- 
ased qubits, e = 0, Hamiltonian evolution via Eqs. Ij28|) 
and (|30j) also neither mixes the spin-0 and spin-1 sub- 
spaces nor the three subspaces of the purely quadratic 
regime. Gradual collapse to one of these subspaces oc- 
curs during the qubit evolution as numerical simulation 
shows. In the linear and non-linear cases there are two 
collapse scenarios: collapse to the spin-0 or to the spin- 
1 subspace, so that the operator (o^ 1 ) + o^ 2 )) 2 is be- 
ing effectively measured. In the purely quadratic case, 
there are three collapse scenarios; we can say that we 
are effectively measuring the operator 
whose eigenspaces are |l) s , |2) B , and the "+" subspace 
{|3) B , |4) B } with eigenvalues —1, 1 and respectively. 
One can confirm using T221 that the measurement sub- 
spaces associated with these collapse scenarios are attrac- 
tive. In other words, small deviations from a collapsed 
state tend toward zero under the evolution 122|) . Simi- 
lar attraction was obtained previously in an analysis of a 
sequence of strong measurements^ 

We characterize below the collapse scenarios by calcu- 
lating the probability of collapse, typical collapse time, 
and power spectrum of the detector current. 

C. Collapse scenarios when e = and U — 
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FIG. 3: Two numerical realizations of two-qubit evolution 
starting from the fully mixed state for a non-linear measure- 
ment, S = —1, by an ideal (n — 1) detector, (a) Evolution of 
Pii . The upper (solid) line illustrates the scenario of collapse 
into the entangled Bell state |l) s , while the lower (dotted) 
line shows a collapse into the orthogonal subspace. The two 
insets show the corresponding spectral densities Si(u>) of the 
detector current, (b) The detector current / averaged over 
the whole time starting from t = is shown for the same col- 
lapse scenarios as in (a). Dashed horizontal lines correspond 
to the normalized averaged currents I 3po , I spi - 



1. Linear and non-linear detection 

Examination of ensemble averaged equations (master 
equations) provides insight into the collapse scenarios. 
These equations can be obtained from ll'L'l) by simply 
eliminating the noise term l|24l) . For both linear and non- 
linear detection, the master equations have the following 
stationary solution 

Pij, st = °> for * + 3 

Pll,st =Pu(0), 

= Pket = PU,st = I [1 - Pfl (0)] ■ (31) 

The solution (pfi(t)) of the ensemble averaged equa- 
tions is a constant in time, (pfi(t)) = pfi(0). Tak- 
ing into account that after some measurement time^ 
t S> T meas all members of the ensemble eventually col- 
lapse to pfi(t) = 1 or px\{t) = it must be that the 



fraction of members that collapse to pf x (t) = 1 is equal 
to (pfi(t)). Thus, one expects that the probability of 
collapse to the spin-0 subspace is /^(O), while the prob- 
ability of collapse to the spin-1 subspace is 1 — pf x (0). 

Extensive Monte Carlo simulations of (1221) were per- 
formed for weak coupling (C between 1/4 and 1). As 
illustrated in Fig. [31 they show that any initial state ei- 
ther collapses eventually into the spin-0 state or collapses 
to the orthogonal, spin-1 subspace within which it per- 
forms oscillations. Within a level of accuracy consistent 
with the number of trials, the probability of evolving into 
the Bell state |1) B coincides with the initial value pf^O). 
In particular, starting from a completely mixed state the 
probability is Pn(0) = 1/4. 

We can monitor the process of collapse of py(t) 
through the Bayesian equation JHJ for a given measure- 
ment result I(t). Collapse is also revealed in the appear- 
ance of different power spectra for I(t) depending upon 
the final state (insets in Fig. E{a)). In the case of col- 
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FIG. 5: Linear weak ideal measurement. State purification 
in case of oscillating spin-1 collapse scenario (pg -> 0, thick 
dashed line) starting from a fully mixed state. The qubits 
signal 2 = 2Rep^, Eq. 125H establishes quantum oscillations 
at Rabi frequency (thick solid line) with slowly fluctuating 
amplitude and phase. Also shown are fluctuations of pf 2 (thin 
dotted line). 



FIG. 4: The two-qubit power spectrum in the oscillating 
scenario in the linear (a) and purely quadratic (b) detection 
regimes. Analytical calculation (dashed line), given by Eqs. 
1471 1 . 1491 1 practically coincides with numerics. 



lapse to the entangled state |1) B , the spectrum is flat 
and coincides with the shot noise spectral density So- In 
the case of collapse to the spin-1 subspace, the qubit per- 
forms oscillations within the subspace. In the linear case 
the oscillations lead to a single peak in the power spec- 
trum (Fig.QJa)) at the Rabi frequency f2 with peak height 
(32/3)S'o and width T = (A/) 2 /4So that is confirmed an- 
alytically below. In the general non-linear case, one sees 
peaks in the power spectrum at the Rabi frequency f2 and 
also at zero and double the Rabi frequency (Fig. Ofa), 
lower inset). 

The average detector current using (0) is (I(t)) = 
J2k h(pkk(t)) t - After collapse, the time average (pkk(t))t 
can be calculated using the stationary ensemble averaged 
solution H31[) - the initial transient period of the order of 
the measurement time can be neglected. In the spin-0 
scenario one finds (p 22 (t)) t = (P33(t))t = 1/2, implying 
(I(t)) = I spa = {h + ^3)/2 = hz- In the spin-1 scenario 
(pu(t)) t = (p4A(t)) t = 1/3 and (p 22 (t)) t = (P3 3 (t))t = 
1/6, so that (I(t)) = I spi = (ii + I 23 + J 4 )/3. In the 
linear case (Fig.|21), the average current in the detector is 
the same in both scenarios ((I(t)} = I spo = I spi = I23) 
since I23 — (Ii + Ii)/2. Therefore the spin-0 and spin-1 
subspaces arc indistinguishable by average currents and 
we have to measure the power spectrum in order to dis- 
tinguish them. In Appendix B we have evaluated the 
typical time needed to accumulate the spectral peak; it 
is on the order of l/r?y 2 . 

In the non-linear case, the final state can be identified 
by the average current in the detector: I spo — 7 23 in the 
spin-0 scenario and I spi = (-Zi+/ 2 3 + /4)/3 in the oscillat- 
ing spin-1 case (Fig.^b)). Thus one expects the current 



signal averaged over sufficiently long times At > 
/ = (1/At) J At I(t')dt' , to be Gaussian distributed with 
average either I23 or {I\ + I23 + ^4)/3 and with variance 
So/2At. We have confirmed this conjecture by numeri- 
cal simulations of the detector current. As a result, the 
typical time necessary to distinguish between the two sce- 
narios with a signal-to-noise ratio of 1 is on the order— 
of r meas = 2S /(I SPO - I SP1 ? = 185 /(2/ 23 - h - hf- 
This estimate is self-consistent since we are working in 
the weak coupling regime (r <§C tt), so that the mea- 
surement time T meas 3> 1 /O. Hence, for At ~ T meas , the 
current / in the oscillating scenario is effectively averaged 
to the value I spi over many periods of oscillation. 

Both in linear and non-linear cases when collapse to 
|1) B occurs, purification of the state is clearly inevitable 
even for a non-ideal detector with r\ < 1. Simulation (see 
Fig. |SJ shows that for an ideal detector, the qubit state 
purifies also in the case of collapse to the spin-1 sub- 
space. (Such purification during oscillation in the weak 
measurement regime is similar to the purification in the 
one-qubit case;** for analogous purification in the case of 
measurement of an oscillator, see Ref. l4l|) . After initial 
relaxation, the surviving nondiagonal matrix elements in 
the Bell basis are found numerically to satisfy 

(Rep^f = pgpu, (Im/4) 2 = P22P33, 

(Impf 4 ) 2 = pf 3 pf 4 (32) 
Impf 4 = Rcpf 3 = Rcpf 4 = (33) 

which imply that the pure-state condition p 2 = p is sat- 
isfied by the matrix p. 

In the weak coupling case C < 1 for a linear ideal de- 
tector (77 = 1) it is possible to characterize further the 
evolution of the density matrix. Since purification of the 
state leads to Eqs. l|3^|) and (|3^|l . of the 8 real degrees 
of freedom in the density matrix of the collapsed system, 
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only 2 remain independent. The system oscillates with 
slowly fluctuating amplitude and phase (see Fig. 0. Tak- 
ing into account the free Hamiltonian evolution, Eq. I)28[l. 
we can write 

z = 2Rcpf 4 = A(t) cos(fif + ip{t)), (34) 
y = 2Impi 3 =A(t)sm(Qt + (p(t)), (35) 
zi = pf 3 - p{\ = -B(t) cos(2Qf + 2<^(t)), (36) 
yi = 2I m( of 4 = -Bit) sin(2ffi + 2<p(t)). (37) 

Here, A(t) = 2^(1-^) and B(t) = 1 - p|, so 
it is natural to regard an d fif) as the two inde- 

pendent variables. The quantities A(t), B(t), and (p(t) 
are constants in the absence of measurement while in the 
presence of measurement they slowly fluctuate in time. 
(The stationary distribution of p\\ is calculated in Ap- 
pendix A.) Eqs. (|3*ljl - i|!r7|) complete the description of 
the collapse subspaces in the case of weak linear ideal 
measurements. 



2. Purely quadratic detection 

In the case of purely quadratic measurement (I\ = 
I4 = /14, 5 = —2), the stationary solution of the master 
equation is 

Pij.st = 0, for i ^ j 

Pn, st =Pu(0), P? 2 ,,t=Pi(0) 

Pis, st = Pu,st = ^ [1 - Pu (0) - P22 (0)] • (38) 

This suggests that with probability pf^O) the system 
collapses under measurement to with probability 

p^ 2 (0) the system collapses to \2) B , and with probability 
1 — /ofi(0) — /*f 2 (0) the system collapses to the remaining 
subspace. 

Numerical simulations verify these three collapse sce- 
narios. The average currents are I23, I14, and (I23 + 
7i4)/2 for the three scenarios so that the current spacing 
is AJ/2. In the case of collapse to the two-dimensional 
{|3) B , |4) B } subspace, the two-qubit system performs 
quantum oscillations as seen in Eq. i|15fl . At times larger 
than the oscillation period (that are relevant in the weak 
measurement regime) we can estimate a typical time nec- 
essary to distinguish between the three scenarios that is 
on the order of T meas w 8S'o/(A/) 2 . 

We checked numerically that the probabilities of the 
three collapse scenarios are given by the initial values 
Ph(0), pf 2 (0)> and 1 - pf x (0) - pf 2 (0) respectively. In 
the case of collapse to the two-dimensional subspace, the 
qubits oscillate, producing one peak in the power spec- 
trum at double the Rabi frequency (Fig.Qfb)), with peak 
height 4S* and width T = (A/) 2 /45 , . 

It is interesting to note the difference between the weak 
measurement results presented here and the hypotheti- 
cal (experimentally difficult) strong measurement case. 



A stro ng, purely quadratic measurement was suggested 
in Ref. |38| as a means of measuring the relative state of 
two qubits for use in error correction procedures. In the 
strong measurement case the qubit Hamiltonian dynam- 
ics is much slower than the measurement dynamics - the 
measurement occurs within a time much shorter than the 
qubit oscillation period r meas ~ Sq/(AI) 2 <C 27r/f2. As a 
result, the two-qubit system collapses to subspaces with 
average current /14 or I23, corresponding to "parallel" 
or "anti-parallel" qubits (Fig.[2J), given in the Bell basis 
by {|2) B ,|4) B } and {|1) S ,|3) B }. Therefore, effectively, 
the operator ai^ai 2 ^ is measured. In the weak measure- 
ment situation considered in this paper, the measurement 
time is much longer than the qubit oscillation period 
Tmeas S> 2tt/Q. During the measurement, a well-defined 
averaged current (^23 + ^i4)/2 appears. Therefore, one 
can distinguish three subspaces leading to effective mea- 
surement of the operator 0^0^ . Contrary to 
naive expectation, weak measurements taken for a long 
period of time are not equivalent to a strong measure- 
ment because of non-trivial interplay between the qubit 
Hamiltonian dynamics and dynamics due to the measure- 
ment process. 

D. Correlator of measurement current and power 
spectrum when e = and (7 = 

In this subsection we derive the current correlation 
function and power spectrum of the detector signal for 
weak measurement of two identical qubits symmetrically 
coupled to the detector. It follows from 1251) and l|26() 
that after collapse to a one-dimensional subspace (to 
in the linear and non- linear detector cases and to |1) B 
or |2) s in the purely quadratic case) the value of z is 
constant, the qubits do not oscillate, and the correla- 
tion function Kj(t) = (I(t)I(t + r)) — (I) 2 contains only 
the noise. The corresponding power spectrum is flat. 
Collapse to a higher-dimensional subspace (the three- 
dimensional spin-1 subspace in the linear detector case 
and the two-dimensional {|3) B , |4) B } subspace in the 
purely quadratic case) leads to oscillations and a richer 
power spectrum that we derive below. 

1. Linear detection 

In the case of collapse to the spin-1 subspace in the 
linear case, we have (I) — I23 and from Eq. H25fl . z — 
pu — P44 = 2Re / o B 4 - Using expression l|26l) . we have 

K^r > 0) = (AI) 2 K z {t) + AIK iz (r). (39) 

First, we calculate the "back-action correlator" K^ z (t) = 
(t;(t)z(t + r)) t for r > 0. In an individual realiza- 
tion, the quantum back action on the qubit state /Oy is 
taken into account explicitly through the Bayesian equa- 
tions ©. In particular, /Oy(t + r) is correlated with 
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noise fluctuations of the current at some previous mo- 
ment t (compare with Ref . l5l|) . Indeed, a noise fluctu- 
ation £(t) leads to a perturbation of the density matrix 
Pij(t) in the immediate vicinity of the time instant t, i.e., 
Pij(t + 0) = Pij(t) + dpij{t). The leading perturbation is 
given by Eq. I^Ijl 

dz(t) ~ -2^ [z 2 (£) - (f&(t)+Pu®)] (40) 
dy(t) ~ -2^ [y(t)z(t) + yi (t)/2] £(t) A (41) 

Do 

where ?/ = 2Impf 3 and yi = 2lmpf i . 

We proceed further by implicitly assuming that the 
long time average (...)( is represented as three sub- 
sequent averages. First, averaging over different real- 
izations of the stochastic process Pij(t') in the interval 
t < t' < t + t is equivalent, assuming stationarity and 
ergodicity, to using the master equations for z(t) and y(t) 



z = -fij/, y = ttz- Ty. 



(42) 



The ensemble dephasing rate is T = r)~ l (A/) 2 /4dV Note 
that the solution to the master equations is linear with 
respect to initial conditions. Explicitly, for z(t) we have: 

z{t) = z(0) exp(-r-r/2)[cos£lT+ (r/2ft)sinft-r] 
-y(0) exp(-rr/2)[fi/0] sin fir 
= z(Q)G z (r)+y(0)G y (T), (43) 

where ft = (ft 2 - T 2 ^) 1 / 2 , and G z (t), G v {t) are the 
corresponding Green's functions. Therefore, the initial 
perturbation that appears due to £(t), propagates to the 
later moment t + r according to Eq. H43|) with "initial" 
conditions dz(t) and dy(£): + r) = dz{t)G z (T) + 

<iy(£) G v (t). The second averaging is the standard aver- 
aging over noise (at instant t) according to (£(£)£(£ = 
Sq/2. The quantum correlation is expressed then as: 

m Z (t+T)) = m5z(t+T)) 

= Al[(pg(t)+pZ(t)~z 2 (t))G z (r) 

+ ( y (t)z(t)+ yi (t)/2)G y (T)}.(U) 

A third averaging over the "initial" condition pfj(t) re- 
mains to be done. 

The calculation for the correlation K z (r) = (z(t)z(t + 
r))t also proceeds in stages, first averaging over different 
realizations for Pij(t') in the interval t < t' < t + r. Then 

(z(t) z(t + r)) - (z 2 (t)) G z (t) - (y(t) z(t)) G y {r). (45) 

In (|39(l the unknown averages (z 2 (t)) and (y(i) z(t)) can- 
cel out, and (y%) = 0, so we obtain 

K t (t) = (AI) 2 (pg+pf 4 )G z (r) = (AI) 2 ^G z (r), (46) 

where we have used stationary values from i|31|) . This is 
valid for arbitrary coupling C and detector ideality rj. 



The corresponding current spectral density Si (to) — 
2 J_ Ki(t) exp(wjr)<iT is readily computed by Fourier 
transform™ to be 



Si(u) = S + 



8 ft 2 (AJ) 2 r 

3 (oj 2 -n 2 ) 2 + rh 



(47) 



In the case of weak coupling, r <§; ft, the spectral peak 
is at ft; it has the same width as in the one-qubit case^S 
and peak height (32/3)ryS'o, which we have confirmed in 
numerical simulations (Fig.Efa)). It is important to em- 
phasize that in an ensemble averaged approach starting 
from a totally mixed state this power spectrum would 
be weighted by 3/4, the probability of collapse into the 
spin-1 subspace. The averaged peak height would there- 
fore be 87760. To derive the higher value of (32/3)7750 one 
must analyze a two-qubit system already collapsed into 
the spin-1 subspace. 



2. Purely quadratic detection 

For purely quadratic detection, the system can col- 
lapse to |1) B or to |2) s , or it can collapse to the sub- 
space {|3) B , |4) B } and perform oscillations. In this last 
scenario, pf x (t) = pf 2 (i) = and (I) = I 23 - AI/2. To 
study the oscillating signal, it is convenient to rewrite Eq. 
(HI in the form I(t) = (I) + {AI/2) zi(i) + £(£) where 
Z\(t) = p§s(t) — p^iif). The Ito equations for the relevant 
density matrix components Zi(t) and yi(t) again decou- 
ple from all other components appearing in Eq. (|2^|l - (|3Tj|l . 
The corresponding master equations read 



ii = -20yi, yi = 2ftzi - Fyi. 



(48) 



Significantly, motion in this two dimensional Hilbcrt 
space is precisely analogous to one-qubit motional, but 
with the oscillation frequency 2ft rather than ft. This is 
natural because quadratic measurement is sensitive only 
to the relative state of the two qubits distinguishing only 
between states with "parallel" and "anti-parallel" spins. 

Applying the same steps as in Eqs. (^UJ-GU; we ob- 
tain the power spectrum 



Si(u>) = S + 



4ft 2 (A/) 2 r 



4ft 2 ) 2 +T 2 u 2 



(49) 



that exhibits a peak at 2ft with height 4t7«S'o and width 
T = (AI) 2 /4r)S just as in the linear case. 30 i 48 i 53 i 54 In an 
ensemble averaged approach, starting from totally mixed 
state the peak height of l|4^|) would be weighted by 1/2, 
the probability that a totally mixed system will collapse 
into the subspace {|3) s , |4) s }. Thus the peak height 
would be 2r]So', it is necessary to analyze a two-qubit sys- 
tem known to have collapsed to the oscillating subspace 
to predict the larger height of 477 So . 
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E. Influence of bias e and two-qubit coupling U 

We now discuss the influence of non-zero energy bias 
e a = £i, = e and two-qubit coupling U for identical 
qubits. Numerical simulations and analytical calcula- 
tions both show that the system still collapses and puri- 
fies. However, differences arise in the collapse scenarios. 

In the case of a linear detector, the system will still 
collapse to the spin-0 subspace |l) s or to the spin-1 sub- 
space {|2) s , |3) B , |4) B }. In fact, the corresponding sta- 
tionary solution, Eq. H31|l . remain unchanged. This is 
because the additional evolution due to e- and [/-terms, 
Eqs. I|29|) and (|30fl . does not mix the two subspaces. 

The characteristics of the spin-0 collapse scenario are 
not affected by the parameters e and U. In particular, 
the resulting power spectrum is still flat with Si (to) equal 
to the detector shot noise Sq. 

In the oscillating, spin-1 scenario the additional dy- 
namics do influence the form of the spectral density 
Si(oj). For [7 = the relevant master equations are 

z = —fly, y = flz — Ty — (e/h)x, 

x = -Yx + (e/h)y (50) 

with x = 2Repf 4 and y, z as in Eq. This is exactly 
the same system of equations as in the one-qubit case 
with a non-zero bias. Thus, the corresponding spectral 
density in the two-qubit case is related to the one-qubit 
spectral density peak S lqb (uj, fl, T, e) = S lqh — S Q , con- 
sidered in Refs. 1301481 

S^^(u i n i r,e) = S + ^§^ b (u,Sl,r,e). (51) 

As in the one-qubit case, finite e leads to a decrease of 
the spectral peak around the Rabi frequency fl and to 
the emergence of an extra peak around zero frequency. 
Note, that the 8/3 enhancement factor is just as in the 
zero bias case, Eq. (14711 . Numerical simulations of Eqs. 
(1231-1(211) confirm this result. 

For e = and U ^ 0, the form of the spectral peak(s) 
in the oscillating scenario can be calculated either by 
numerical simulations using Bayesian Eqs. (|23l) - (|30[l or 
through an ensemble average method developed in Refs. 
I35H31 Increasing U from zero first suppresses and then 
splits the spectral peak at the Rabi frequency fl. 

In the case of purely quadratic detection, qubit biasing 
so that £ ^ while U — will lead to mixing between 
the \2} B and {|3) B , |4) B } subspaces as seen from Eq. (|2"31) . 
For small bias e this will take the form of rare switching 
events whose rate and influence on the spectrum will be 
discussed in the next section. 

If U 7^ while e = 0, one notices that the density ma- 
trix components Zi(t) = p B 3 (t) — pf^t), yi{t) = 2Imp B 4 , 
and x = 2Kepf 4 obey equations analogous to (|50|l in 
which fl is replaced by 2fl, z by Zi, y by y\, and e by 
2U. Thus, the form of the spectrum is described with 
the same functional form, as in Eq. (|51|l 

S 2 I qb - quadr {u J ) = So +~S) q \uj,2 fl,T, 2 U). (52) 



The result is a decrease of the spectral peak at 2fl seen in 
Eq. O and the appearance of a peak at zero frequency. 
This finding is confirmed by numerical simulations^ 



IV. QUANTUM ZENO-LIKE STABILIZATION 

The quantum Zeno effect is usually understood to arise 
in the case of projective measurements taken more and 
more frequently. In our case of continuous measurement 
the detector is always coupled to the system, so the ap- 
proach to the quantum Zeno regime corresponds instead 
to the limit of stronger and stronger coupling. Since 
it is impractical to realize strong coupling detectors in 
many systems (the dimensionless coupling to each qubit, 
C ~ T/fl, is typically small), one might not expect quan- 
tum Zeno physics to arise. However, in our system of 
two identical qubits, a small perturbative asymmetry be- 
tween the qubits, such as when Afl = fl a — fib =/= 0, can 
generate slow transitions between states that were uncou- 
pled. The effective coupling for this dynamical process, 
C e ff ~ r/AQ, can be very large, leading to a Zeno effect. 

From this perspective, the symmetric-qubit situation, 
AO = 0, constitutes the Zeno limit of infinitely large cou- 
pling C e ff, the transition becomes completely inactive, 
and the quantum state becomes "frozen" . The asym- 
metric situation with non-zero Afl implies finite effective 
coupling Ceff- The non-zero Afl makes the quantum 
state to "decay" from one subspace to another. But, 
since the coupling is large, the decay of the quantum state 
looks like an abrupt switching event, with duration on the 
order of the measurement time T meas ~ 1/r, rather than 
gradual evolution of frequency Afl. The switching event 
can be identified experimentally by measuring the power 
spectrum, for instance. From an ensemble averaged point 
of view, at long times the decay becomes exponential with 
a rate that is much smaller than T. 

Numerical simulations using Bayesian equations © 
confirm that the two-qubit density matrix p experiences 
rare, abrupt switching between different subspaces when 
the qubits have a slight asymmetry like a non-zero Afl. 
These switchings are apparent in Figs. and [7] Switch- 
ing also happens for identical qubits when there is small 
difference of couplings C a , to the detector (due to dif- 
ferent values for A/ a and AIj,) or different coupling to 
an environment. 

For linear detection the switching occurs between the 
spin-0 and spin-1 subspaces while for quadratic measure- 
ment we will have three-terminal switching among the 
three collapse subspaces described above. The switching 
events can be directly observed through measurement of 
the intermediate time scale power spectrum and/or av- 
erage detector currents. (Such intermediate time observ- 
ables are meaningful if the measurement time needed to 
accumulate the spectrum is much smaller than the "life- 
time" of the subspaces.) If measurement is performed 
for a sufficiently long time, it will lead to averaging of 
the spectrum and may also lead to the appearance of an 
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FIG. 6: Linear detection. Switching between spin-0 and spin- 
1 subspaces as seen in the time evolution of pfi(t), in case of 
different Rabi frequencies. The parameters are: C = 1/4, and 
Afi/fl ~ 2.2 x 10 -3 . The inset shows numerical histograms 
(in logarithmic scale) of switching time distributions in agree- 
ment with the asymptotic switching rates Eqs. 158H . 163H . 



additional telegraph noise peak at zero frequency anal- 
ogous to the regime of strong coupling in the one-qubit 



1. Linear detection 

For linear detection (S = 0), assuming weak coupling 
An <C r <C n and unbiased, non-interacting qubits 
£a = £b = U = 0, a system of equations for the quantities 
Pli - P22, 2Impf 2 , pf 2 - p B 4 , pf 3 - pf 4 , and 2Impf 4 can 
be extracted from Eqs. J22J|-|2II) an d To change 

them into a system of master equations for ensemble av- 
erages, we simply set the noise term to zero £ = in 123( 1 . 
Starting from the entangled initial condition p B i(0) = 1 
we find pf x (i) = 1 - J* dt' 2 Impf 2 (*')- Formal inte- 
gration of the equation for 2Impf 2 (see Eqs. I)23fl. JHSJ) 
leads to the exact relation 



^2 r t 
2' ,„ 



dt"e rt " (pfUf") -*&(*"))• (55) 



Solving our system of master equations to zeroth order in 
An and substituting into l|55|l . we obtain an approximate 
solution for the "survival probability," pfi(t): 



case 



30,31 



pfl(t) = l + 



(An) 2 (An) 2 (An) 2 

2r 2 2T~ 2T 2 



exp(-rt). (56) 



A. Perturbation to the qubits 

To find the switching rates analytically, we have used 
the master equation, which reads of Eqs. I(23|l (|29ll 



dpfj 






dpfj 




dp% 




dpfj 


dt 


dt 


^ 

meets 


dt 


H 

n 


dt 


€ 


dt 



with the noise term I24H omitted from dp B J dt\ . It 

1 1 V <meqs 

now acquires additional terms. When the frequencies of 
the two qubits are slightly different while coupling to the 
detector is symmetric, I 2 = ^3, we have in addition to 
Eq. |2EJ) 
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and in addition to Eq. (|2"§|) 
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For times £ <C 1/r, one can expand the exponential 
and find that there is no linear term in the i-expansion 
of 1561) . This fact makes quantum Zcno effect physics 
possible The first few terms in the small i-expansion 
of H56[) actually coincide with the exact small i-expansion 
of pfi(t) derived directly from the master equations by 
differentiation, using the initial conditions 



P?i(t) = 1 



(An) 2 



t 2 



(MTfr 
12 



(57) 



-(An) 2 /4 



The ^-coefficient turns out to be 
— Tr[/5 2 7ig B ] + TtIpHqbpHqb], i-e- it is determined by 
coherent (Hamiltonian) evolution alone consistent with 
the discussions in Refs. I ll5l l(i One might speculate that 
the system has some memory at the time scale t <C 1/r, 
"forgetting" its history once t > when the exponen- 
tial term in Eq. (|56|l drops out. At still longer times, the 
"decay" from the Bell state to the oscillating subspace 
eventually becomes exponential with rate 



B^O — 



(An) 2 
2T 



(58) 



as demonstrated by switching time distributions shown at 
the inset of Fig-EI Note, that the switching rate in l|58|) 
is similar to the one-qubit case^S with the replacement 

n -> An. 

An analogous derivation for the case of different qubit 
biases e a can be performed. We set e = (e a +e b )/2 = 
since effects of having e 7^ (while e a = £&) have 
been considered above, in Sec lIIIEl We stipulate that 
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\Ae\/h <C ft but permit |Ae|/fi, and T to be of the same 
order. From (J2SJ, (|54l) we find 

pfi(t) - -(Ae/ft) Impf 3 - ii 3 (t)/2, (59) 

where 2:13 (t) = pfi(t) — pf^it). The second equality in 
()59J) was obtained by solving a system of nine equations 
extracted from Eqs. 12311 . 128fl . and l|54|) . using the initial 
condition pfi(0) = 1. These nine equations quickly re- 
duce to a system of three equations for 2 lmpf 3 , 2 Repf 4 , 
and X13 = jO^— p^j, which can be solved perturbatively in 
Ae/h. Integrating JSHJ, we obtain pf x (<) = (l+xi 3 (t))/2. 
An approximate solution for p B x (t) is 



1 — Dr—t ~\ 

£ 2 2 



exp(-R/2)x 



c e cos ft?; — [2 — c £ ] sin fit 

2ft 



, (60) 



where b F 



r 2 /[ft 2 



{Ae/fif/p 2 + (Ae/ft) 2 ], c £ 
(Ae/Ti) 2 ] - 1, and ft = + (Ae/?i) 2 - T 2 /4. As in 

(|57|l . there is no linear term; the evolution is quadratic in 
t at small times permitting quantum Zeno effect physics^ 
The decay is exponential on a coarse time scale t > 2r _1 
with switching rate 



B^O 



(Ae/h) 2 r 
2ft 2 ' 



(61) 



Numerical solution of the master equations for the den- 
sity matrix confirms these results. 

We have computed the rate of switching from the spin- 
Bell state |1) B to the spin-1 subspace. To compute the 
rate of reverse switching, we notice that the stationary 
solution of the master equations (Eqs. (|2*3^l - l(3*n|l . (p^jl . 
and H54f) with the noise £ set to zero) has the form 

Pht = Hi,* = %/ 4 - ( 62 ) 

This form is completely mixed because (|53f) and l|54|) 
cause transitions between the spin-0 and spin-1 sub- 
spaces. One infers that the system should spend on av- 
erage 1/4 of the time in the state |1) B , suggesting the 
stationary state condition (1/4)Tb^o — (3/4)ro^s = 



r 



O^B 



B^O 



/3. 



(63) 



The numerical histograms of switching time distributions 
(see inset of Fig. EJ), which were accumulated in a Monte- 
Carlo computation, confirm the long-time exponential 
decay of the "survival" probability pfi(t) as well as the 
asymptotic rates l{55j) . l{oTl) . and l(6"3*|) . 

If a state's lifetime in the oscillatory subspace, to = 
1/To^b, is sufficiently long, it is possible to measure the 
spectrum described by Eq. (|47|l . The necessary condition 
for this is to ^ T meas where r meas ~ 1/r is a typical 
measurement time needed to achieve signal-to-noise ratio 
close to unityi2& For small efficiency ij we estimate (see 
Appendix B) 



(64) 



32r?7 2 



while for 77 close to unity the measurement time may be- 
come an order of magnitude larger than l|64|l . Physically, 
the lifetime will satisfy to T meas if one can tune close 
to the symmetric point of identical qubits (small Aft) 
and if the efficiency 77 is not too small. 

The spectrum accumulated for times much greater 
than the lifetime to gets averaged by switching events. 
Since the system resides in the spin-0 subspace 1/4 of the 
time and in the spin-1 subspace the remaining 3/4, the 
spectrum will take the form l|47|l with a factor of 2 instead 
of 8/3 (i.e. the peak-to-pedestal ratio will be 877). There 
will be no additional peak at zero frequency because the 
averaged currents in the two scenarios coincide. 



2. Purely quadratic detection 

For a purely quadratic detector (S = —2) in the case of 
r 3> Aft 7^ 0, Eq. (|53|l produces switching transitions be- 
tween states |1) B and \2) B while the oscillating state in- 
side {|3) B , |4) B } does not switch. The stationary solution 
to the master equation is similar to (|38|) , but now p B x st = 

P22 st — 5 iPii (0) + P22 (0)] ■ The new measurement sub- 
spaces are the minus subspace "— "= {|l) s ,|2) s } and 
the plus subspace "+"= {|3) B , |4) B }, where the names 
are based upon l|10|) - (|13fl . The switching rates are given 
by 



Tli) 



\2)i 



= r.. 



(Aft) 2 

2r 



(65) 



derived by solving master equations to lowest order in 
Aft as in the derivation of Eq. I|58l) . Switching between 
|1) B and \2) B , and therefore between currents I23 and 
7i4 in Fig.[3 produces telegraph noiseS^ As a result, if 
we measure long enough (t ^> I/T^b^^b), the power 
spectrum will exhibit a peak at zero frequency 



Si(uj) = S + 



1 



(Ai) 2 r 

(Aft) 2 1 + [^r/(Aft) 2 ] 2 



(66) 



in full analogy with the one-qubit case. 

Actually, it is possible to derive the power spectrum 
for arbitrarily large Aft: 



St(u>) = S + 



(Aft) 2 (A/) 2 r 



(Aft) : 



T 2 L 



(67) 



which reduces to (|66|l for small Aft. Indeed, after col- 
lapse to the "-" or "+" subspaces, Eqs. {T1J, CSJ, 
l)23[l . and l|24|l produce weakly perturbed quantum os- 
cillations within the subspace and no switching to the 
other subspace. The quantum oscillations in the mi- 
nus and plus subspaces correspond directly to one-qubit 
dynamics with oscillation frequency Aft and frequency 
2ft = (ft a + ft;,) respectively. After a collapse in the "— " 
subspace the spectrum will be given by (|6Tj) . If collapse 
has happened in the "+" subspace the spectrum will be 
given by Eq. I|49|l with a peak at 2ft. In both cases the 
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FIG. 8: Current levels for asymmetrically coupled qubits, 
AI a 7^ Linear and quadratic response. 
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FIG. 7: Quadratic detection. Switching between Bell sub- 
spaces |2) s and {|3) s , |4) s } of Eqs. jTUt - ifT-fll as seen in 
the time evolution of p^C*): due to a small qubit bias. |l) s 
is unaffected as seen from 1291 . The parameters are: C = 1, 
and e/fi ~ 2.5 x 10" 2 . 



maximum peak-to-pedestal ratio at a non-zero frequency 

iS 4:7]. 

Measuring the spectrum on a sufficiently long time 
scale will constitute an (effective) measurement of the 

operator a^a^ whose eigenspaces are the "— " and "+" 
subspaces. 



As we mentioned in Sec lIII El if instead of an asym- 
metry in frequency, the qubits are identical and biased 
£„=£(, = £ C T C 11, under quadratic measure- 
ment switching will occur between the \2) B state and the 
{|3) B , |4) s } subspace (see Eq. (J2HI))- The spin-0 subspace 
|1) B will be unaffected as shown in Fig.[7| 

The asymptotic switching rates (relevant at times t > 
1/F) are derived by study of master equations. The sta- 
tionary solutions are given by (|31|l and the switching rate 
by 



2e 2 

r|2)«^{|3)«,|4)B} = "^2" r - 



(68) 



= 



Since pf 2 , s t — P33,st — Pii,st — 1/3, Pii, s t — Pf^j,st 
is a stationary solution of the master equation within 
the spin-1 subspace, the stationary condition must be 

( 1 /3)r| 2 )-B- t {|3)- B ,|4)- B } - ( 2 / 3 ) r {|3> B ,|4) fl }->|2) B = 0, im- 
plying 



{|3)B,|4>S}->|2>' 



1 



*{|3> B ,|4} S }- 



(69) 



When the system collapses to |1) B , the power spectrum 
remains flat Si (to) = So- In the case of switching between 
1 2) B and {|3) s , |4) B }, we can average the spectra as in 
the linear case above and add a telegraph noise peak 
caused by switching between current levels Iu and (J23 + 
ii 4 )/2. The result is 



Si{u) = S a + - 



4n 2 (AI) 2 T 



Note that the peak-to-pedestal ratio at 2SI is multiplied 
by a factor 2/3 as a result of the averaging. This re- 
sult can be reproduced also in the ensemble averaged 
approachi^2& 

To conclude this section we mention that a non-zero 
qubit interaction U does not mix the measurement sub- 
spaces, neither in the linear nor in the quadratic case. 
Correspondingly, calculations show that neither the sta- 
tionary solutions of the master equations nor the switch- 
ing rates derived above will be affected. Finite U will 
only affect the form of the spectrum in the oscillating 
scenario as described in Sec lIII El 



B. Slightly different couplings 

For identical qubits with slightly different coupling to 
the detector, I 2 ^ I3 as indicated in Fig.|SJ This asymme- 
try introduces evolution between collapse subspaces. The 
presence of the detector provides quantum Zeno stabiliza- 
tion, and so instead of gradual evolution, there is rare, 
abrupt switching between the subspaces. It is convenient 
to introduce the current asymmetry dil = (I 2 — l3)/2 = 
(AI a — AIb)/2 and to define the current difference as 
AI = AI a - 5x1 = AI b + hi = (A/ a + AI b )/2. The 
current asymmetry parameter is 8\ = 5\I / AI. Defining 
the difference in coupling by AC = C a — C and the aver- 
age coupling by C = (C a + C )/2, we note that in cases of 
small coupling asymmetry, 5i 



AC 
4C ■ 



1. Linear detection 

If 61 ^ 0, the following additional terms arise in the 
master equation in the case of linear detection 
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If the system starts in the entangled initial condition, 
pfi = 1, the exact solution to the master equation for 
small times is 



pfx(i) - l-2(TS 2 )t + A{T6 2 ) 2 t 2 



(72) 



The solution contains a linear term setting a time scale 
(TS 2 ) _1 > r _1 . For p B (*) close to unity, the asymptotic 
rate of switching from the spin-0 Bell state |1) B to the 
spin-1 subspace is found to be 



B^O 



2T8l 



r]8S 



AC 
C 



= o hr r, (73) 



where the qubit detector couplings are slightly differ- 
ent, AI a 7^ Alb- The reverse switching rate is given 
by Eq. 

To understand Eq. 1)73(1 ■ we note that measurement for 
a short interval 8t cannot distinguish between the current 
levels h and I3 when {h^h) 2 < £>/ where Dj = S /2<Si 
is the variance of the measured noisy current. For a sys- 
tem initially in the spin-0 state 1 1) B ~ |2) — 13) weak mea- 
surement tends to collapse the system towards either |2) 
or 1 3), thereby diminishing the amount of entangled state 
|1) B . Standard decoherence^ in the measurement basis 
P23 — ► /02 3 exp(— {1 2 — h) 2 /%D[) leads in the Bell basis 
to pf^t) ~ 1- [(J 2 - /a) 2 /^] <5t in agreement with 
Note that the rate 173|) is one half of the average classical 
information acquisition rate 39,47 (r c ; = (7 2 — I3) 2 / 't?4iSo) 
for a weakly measured hctitious two level system, as ex- 
pected. 



2. Quadratic detection 

In the case of quadratic detection, if I2 and I3 are 
slightly different while I\ remains equal to J4, then the 
two qubits have slightly different couplings and also the 
detector deviates slightly from being purely quadratic. 
The additional terms in the master equations take the 
form (|71|l but with the index "2" interchanged with the 
index "4" everywhere in the first matrix. As a result, the 
state 1 2) B is unaffected, while rare switching develops 
between the state |1) B and oscillations in the subspace 
{|3) B ,|4) B }. The small time expansion of pfi(t) takes 
the same form l|72|l and thus the rate ^\i) B —>-{\3) B ,\4:) B } is 
the same as in the linear case, Eq. (|73H : 



is r. 



In the case of 



{|3)*,|4)*}-|i> fl - 5 r |i> B ^{|3> B ,|4)B}- 
collapse to |2) B , the spectrum remains flat, 5j(w) = So- 
Otherwise, for measurements over a time interval too 
short for switching, the system will reside in |1) B leading 
to Si(ui) = S or in the subspace {|3) B , |4) B }, leading to 
Eq. Q49J1 . For longer times, the spectrum will be the av- 
erage of Si(ui) = So and the spectrum given by Eq. l|4*§|) 
together with a telegraph noise peak: 
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(75) 



The peak-to-pedestal ratio at 2ft will be (8/3)77 rather 
than Arj as a result of the averaging. 

If the detector has I\ ^ I4 while I2 — I3, switching 
occurs between the state \2} B and the oscillating sub- 
space {|3) B ,|4) B }, while state |1) B is unaffected (as in 
FigEJ . This slightly non-quadratic mode of detection is 
"orthogonal" to the I2 7^ I3 mode in the sense that the 
switching rate would not depend on a small additional 
I2 7^ ^3- asymmetry. The rates are 
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r {|3) B ,|4)B}^|2)f - 2 r |2) B ^{|3) B ,|4)B}- ( 7 7) 

The power spectrum for large times will be given by Eq. 
(17511 with the substitution (I2 — I3) — > (Ii — 14). 

When I2 =/= I3 and I\ 7^/4, we combine the two cases 
h 7^ ^3, h = I4 and I 2 = I3, h 7^ J4. Based on 
the analysis above, there should be switching |1) B 
{|3) B , |4) B } <-> |2) B but no direct switching between |1 
and |2) B . The stationary solution of the master equation 
i s Pu,st = 1/4, pf^j.st = 0, so only one half of the time is 
spent in the oscillatory state, and therefore the peak-to- 
pedestal ratio at 2ft is 2ry rather than Arj. The current 
spectral density in this case, including a low-frequency 
peak due to telegraph noise in a 3-terminal system, is 
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where Ti = r [i)f ->{|3> 
are given by (|74|) and 



and r 2 



|2)B-{|3>* |4)*} 



r |l) B ^{|3) B ,|4)B} 
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(74) 



Since pf l st = pg st = p^ st = 1/3, p B ^ j st = is a 
stationary solution to the master equation within the 
subspace {|1) B , |3) B , |4) B }), the reverse switching rate 



C. Environmental dephasing 

We have studied the effect of environmental dephas- 
ing, modeling it with two small dephasing rates 7 a and 
7b acting separately on the qubits. This adds an extra 
dephasing contribution 7^ to 7^ in Bayesian equation 
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We assume that the coupling to the environment is much 
smaller than the coupling to the detector, 7 a ,7h <C T. 
Numerical simulation again shows that abrupt switching 
events take place. However, now the density matrix p 
remains slightly mixed even for an ideal detector. 

In the linear case we have switching between spin-0 
and spin-1 subspaces. Analyzing the resulting master 
equations, we derive in the case of a linear detector the 
switching rates 



p7<j,7f> 
1 B^O 



qp7a,7f> 

01 o->s 



( 7a + 7b )/2. 



(80) 



In the quadratic case the switching involves transi- 
tions |1) B <-> {|3) s , |4) B } <-> \2} B , with again no di- 
rect switching between |1) B and \2) B . The rates are 



\l) B ^{\3)' 



,|4}«} - r | 2>i 



>{|3) fl ,|4> B } = (7a+7b)/2, and 
the reverse switching rates follow from the stationary 
state. 



V. CONCLUSION 

In this paper we have investigated in detail the non- 
trivial interplay between the internal two-qubit Hamilto- 
nian evolution and two-qubit dynamics under weak con- 
tinuous measurement described by a Bayesian formalism. 
The gradual collapse of the system's quantum state takes 
a finite measurement time T meas ~ On a time scale 
greater than 1/T, the measurement effectively looks like 
orthodox projective measurement in the Bell basis. A 
linear or non-linear detector measures the operator of 
total "spin", (ct >(1) + ~^ { - 2) ) 2 . The system's possible col- 
lapse subspaces are the spin-0 Bell state, |1) B and the 
spin-1 subspace within which the system executes oscil- 
lations. For purely quadratic measurement, the operator 
a y a v + Oz Oz is measured so that the collapse sub- 
spaces are the spin-0 state, the other Bell state, \2) B , and 
the oscillation subspace {|3) B , |4) B }. 

Measurement leads to spontaneous entanglement 
which can be identified by a distinct power spectrum 
and/or distinct value of the average current of the detec- 
tor. We obtained analytically the power spectrum of the 
current, Si(u>), in different situations. In particular, in 
the case of a linear (purely quadratic) detector response 
for non-interacting identical qubits the spectrum exhibits 
a Lorentzian single peak at frequency Q (20) with peak- 
to-pedestal ratio of (32/3)r/ {Arj) respectively. We also 
obtained analytic expressions in the case of qubit bias- 
ing and in the presence of qubit-qubit interaction. For a 
general non-linear measurement the spectrum can be cal- 
culated using numerical simulations of the measurement 
process. 



Various small imperfections of the two-qubit system 
such as asymmetry of the Hamiltonian and/or asymme- 
try in the coupling to detector cause transitions between 
collapsed states. Because the coupling to the detector, 
although weak compared to most of the energy scales 
in the Hamiltonian, can be strong compared to the small 
imperfection, a regime of quantum Zeno-like stabilization 
arises. The system does not drift gradually from one col- 
lapse subspace to another; instead it remains stabilized 
in one subspace and experiences rare, abrupt switching 
events to other subspaces. For very long observational 
times, because of the switching, one observes only aver- 
aged dynamics in the power spectrum. 

The authors gratefully acknowledge the support of the 
Packard foundation (R.R. and A.M.) and of NSA and 
ARDA under ARO grant W911NF-04-1-0204 (A.K.). 



APPENDIX A: CALCULATION OF THE 
STATIONARY DISTRIBUTION FOR pg 

In this Appendix, we use a Fokker-Plank equation to 
calculate the stationary distribution for p B 2 and its mo- 
ments in the case of linear detection in the oscillating 
scenario of Fig.0and Eqs. l(51 )l -l|37 jl . 

An Ito stochastic equation for p 22 can be derived from 
Eqs. (J^HJnEEl, and the parameterizations l|3*I )) -[|3T )l . 
that depends only on the two independent variables, 

x{t) = pf 2 (*) and 

dx A/ 

— = 2 — y/x(\ - x) (1 - 2x) cos(m + <p(t)) f(t) 
at So 

-2T[x- {1 ~x) cos 2 (nt + tp(t))]. (Al) 



In the weak coupling regime one can average over the 
oscillations with frequency O. Then (cos 2 (ftt + <p(t))) — > 
1/2 and 

dx AI 



— = — V2x(l - x) (1 - 2x) fc(f) - T (3a: - 1), (A2) 
at bo 

where, £i(i) is an effective white noise with the same 
spectral density. 

A Fokker-Plank equation for the distribution P(x, t) 
follows^i from Eq. J22J 

^ = -» [B ( a )P( a> *)] + ^^ [C -(,)P(, f *)] 

(A3) 



with B(x) = -r(3a: - 1) and C(x) = %^y / 2x(l - x). A 
stationary solution of Eq. (|A3(I can be derived^ taking 
into account that restricted variable and 

therefore the "probability current flow" is zero. 
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exp 


2(1- 2x) _ 



x 1 / 2 (l-2x) 7 / 2 



(A4) 



with normalization constant Ao = 6yV/e. Eq. (|A4|) 
implies that p B 2 actually fluctuates between and 1/2, 
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which is confirmed by Monte-Carlo simulations, Fig. [3] 

- 1/2 

For the moments, (x n ) = J Q x n P(x) dx, one obtains 



(a;) = 1/3, (x 2 ) = 1/8, 



(A5) 



*yfa x Fi{\-n,-\,-\) (A6) 



8V^T(n+ |) 
15 T(n-2) 



J \*i(3-n^;-±) 



2l/2— n 



where 1 Fi(a, 6;z) is the confluent hypergeometric func- 
tion. 



APPENDIX B: NOISE OF THE NOISE AND THE 
MEASUREMENT TIME TO ACCUMULATE THE 
SPECTRUM 

Here we calculate the measurement time to distinguish 

between a flat spectrum (shot noise) and a spectrum that 

possesses a peak on the top of the shot noise pedestal. To 

make an estimation we first consider a situation with a 

narrow band (A/ -C /, / = lu/2tt) frequency filter. The 

average power for a finite time At (> 1/A/) is a fluctuat- 

— 2 
ing quantity, V — -57 J_ 2 ^t dt (I(t) — Jo) , that depends 

2 

on At and on realizations of the process. Denoting the 
fluctuating part of the current as I(t) = I(t) — I we 
perform averaging over realizations using the well known 
relation to the noise power spectrum (see, e.g. Ref. \E^l : 
{I 2 {t)) r ~ Si(uj)Af. Therefore, averaging of V gives 



At 



dt (I(t) - Io) 2 ) -S/HA/. (Bl) 



The second moment, ("P 2 ) r , can be evaluated for a Gaus- 
sian random process I{t), using stationarity: 



(V 2 



1 

At 



dt (I(t) - I f 



1 r°° 

^AtJ d ^( T ) + (^MA/) 2 , (B2) 

where Kj2 (r) is the correlation function of the "power 
signal" I 2 (t) - (I 2 (t)) r 



K p {t) 



P(t)~(I 2 ) r P(t + r)-(P), 



(B3) 



The variance of the fluctuating power, Eq. (|B2J| . is ex- 
pressed through the "second" spectral density Sj2 (u;) = 



2 J_ QO dr Kp(r)e WT , which for a Gaussian process I(t) 
is related to the ordinary one^i In particular, 5/2(0) = 
2[S7(w)] 2 A/ and we obtain 



AV 2 ^{V 2 ) r -{V) 2 r = — [S I {uj)] 2 Af 



(B4) 



For a more general frequency filter a{f) the power 
is distributed with average and variance respectively 
given as (V) r = J °° a(/)S;(w)d/ and AV 2 = 
M Io° a2 (/) [ S i(u)} 2 df that generalize Eqs. (TrTTi . llmll . 

We are interested in measuring I(t) to distinguish be- 
tween two hypotheses: (1) flat power spectrum Si(lo) = 
So and (2) power spectrum with a peak Si(u>) = So + 
Speak(f)- Assuming Gaussian distribution of the power 
signal we can estimate the measurement time T meas 
at which the variances AV 2 [a(f)], A"Pf [a(/)], become 
small enough to distinguish between the corresponding 
averages (Vi) r , (T-^r- For a "signal-to-noise" ratio of 1, 
Tmeas is expressed as (see also Refl26T) 



y/ AT> 2 [a(f)] + APi[a(f)} 



[fo°° a (f) S peak(f)df]' 



(B5) 



Optimization over the filter requires ST meas /Sa(f) = 0, 
which leads to the filter form 



«(/) 



Speak (,/*) 



l + a{l + S peak (f)/S ) 2 



(B6) 



that is close to the peak form itself while further opti- 
mization on the parameter a is needed. However, for a 
small peak height, h ps S pea k(Q) /So <C 1, the result is 
independent of a and reads 
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Th 2 ' 



(B7) 



Note that small height translates to small detector ideal- 
ity (for one qubit h < At] while for two qubits h < 32/7 / 3) . 
To get an exact result for large or small height, one em- 
ploys integration by residues and find a cumbersome ex- 
pression that can only be managed numerically. Numer- 
ically, the calculation suggests that a op t is of the order 
of 1. An approximation with accuracy better than 30% 
even at large 77 is given by 
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1 



3/i 

T 



77^ 

"64~ 



(B8) 



which means that r meas may become an order of magni- 
tude larger than the simple estimate, Eq. (|B7fl . 
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